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ABSTRACT 



Context. Star-forming small galaxies made out of collisional debris have been found in a variety of merging systems. So far only a 
few of them are known in Ultraluminous Infrared Galaxies (ULIRGs) although they show clear signs of interactions. Whether external 
star formation may take place in such objects in an open question. 

Aims. The aim of this paper is to identify and characterise the physical and kinematic properties of the external star forming regions 
in a sample of ULIRGs, including TDG candidates. The likelihood of survival of these regions as TDGs is also evaluated. 
Methods. The analysis is based on optical Integral Field Spectroscopy (IPS) and high angular resolution HST imaging. 
Results. We have found that the presence of external star-forming regions is common with 12 objects being identified in 5 ULIRGs. 
These regions show a large range of dynamical mass up to 1 x 10'" Mq, with average sizes of ~750 pc. In addition, the line ratios (H ii 
region-like), metallicities (12-1- log (0/H) ~8.6) and Ha equivalent widths (34 - 257 A) are typical of young bursts of star formation 
(age ~ 5 - 8 Myr), and similar to those of the TDG candidates found in less luminous mergers and compact groups of galaxies. The 
extinction corrected Ho- luminosity of these young bursts leads to masses for the young stellar component of ~ 2 x 10* - 7 x 10** Mq. 
The likelihood of survival of these regions as TDGs is discussed based on their structural and kinematic properties. Particularly 
interesting is our finding that most of these systems follow the relation between elfective radius and velocity dispersion found at lower 
(globular clusters) and higher (Elliptical) mass systems, which suggests they are stable against internal motions. The stability against 
forces from the parent galaxy is studied on the bases of several criteria and a comparison of the data with the predictions of dynamical 
evolutionary models is also performed. Five regions out of twelve show High-Medium or High likelihood of survival based on all the 
utilised tracers. Our best candidate, which satisfy all the utilized criteria, is located in the advanced merger IRAS 15250-1-3609 and 
presents a velocity field decoupled from the relatively distant parent galaxy. 

Key words, galaxies: active — galaxies: interactions — galaxies: starburst — galaxies: dwarf 



1. Introduction 

, Ultraluminous Infrared Galaxies (ULIRGs) are defined as ob- 
' jects with an infrared luminosity L(8 - lOOOyum) ^ 1O'^L0 
(e.g. Sanders and Mirabel, 1996). This huge luminosity is be- 
lieved to be caused mainly by star formation, although the ex- 
istence of an AGN cannot be discarded and could be the dom- 
inant source of energy in a sma ll percentage of these systerri s 
(for reviews see ISanders & Mirab el 1996; Lon sdale et alJl2006h . 
ULIRGs are systems rich in gas and dust (e.g. iEvans et al.ll2002L 
lEvans. Surace. & Mazzarellal l2000l) and all of them present 
emiss i on Unes in the optical (e.g. Veilleux. Sanders. & KimI 
Il999t iKim. Veilleux. & SandersI 1 19981) . Observational studies 
about the morphology of these systems indicate that most (if 
not all) of the m show sig ns of m ergers and interaction (e.g . 
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Cui etalJl200'i] 



Veilleux et al J 120021) while theoretical simula- 



tions of mergers indicate that ULIRGs constitute a well defined 
subgroup within mergers in general. That is, in principle only 
major mergers of two gas-rich spiral galaxies are able to ex- 
plain the large luminosity observed in UL IRGs and only dur- 
ing a short period of time (~ 50 Myr, e.g. iMihos & HernquistI 



ll996l:lBekki & Shioval200l . In this context, there are some em- 
pirical evidences suggest ing that ULIRGs a re mergers of two 
or more (sub)-L* spirals dCoUna et"ani200ll), that could be the 
progenitors of interm ediate-mass ellipticals dGenzel et al.ll200il : 
Tac coni et alj 2002h . ULIRGs could also evolve into quasars 



(Sanders et al 



19881) . at least in t hose mergers wher e the pro- 



genitor galaxies are more massive (ICoIinaetaI.ll2001h . 

One issue regarding ULIRGs that deserves attention is the 
presence of knots and condensations of star formation outside 
the nuclei. On the one hand, at the lowest masses and smallest 
sizes ran ge, we find the so-called Super Stellar Cl usters (e.g. 



Scoville e t al. 2000; Shioya et al. 2001; Su race & Sa nders 20001 
I999h . These systems, with typical masses in the range of 
10"* - 10^ Mq and sizes of r ~5 pc, are be lieved to be the 
proge nitors of today globular clusters (e.g. ISchweizer et al.l 
Il996 f) and, in addition to ULIRGs, they have been found in 
a wide range of environmeiits such as starburst galaxie s (see 
for instance iMeurer et alj 1 19951 and iMelo et all l2005h. the 



space between galaxies in compact groups ([Gallagher et al.l 
120011: iMendes de Oliveira et alj |2004|) and speciall y in merg- 
ers l ess luminou s than those studied h ere (see 'Zepf et alj 
19991: IWhitmore et alJ 119991; lAIonso-H errero et al. ' 200l 
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Knierman et al. I l2003r iBastian et alJ 1 120061 and others). On the 
other hand, at masses ~ 10^-10^ Mq we find the so-called 
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Tidal Dwarf Galaxies (TD G). Their existence were already 
proposed in the 50's when IZwickvl suggested that the tidal 
forces in interacting galaxies could create self-gravitating 
objects made up from the debris of the interaction that, per- 
haps, would become small galaxies. From the observational 
point of view, this kind of objects, with masses and sizes 
similar to d warf galaxies and large quantities of gas (e.g. 
iBraine et al.l [2001) have already been foun d in less lumi- 
nous i nteracting or rn e rging galaxies ( Due & Mirabelj |1994 
1998[ iDuc et al.' '2000^, 'Hibbard et al.' '2001'; 'WeilbachereLalJ 
200d IWeilbacher, Due, & Fritze-v. Alvensleben 2003^ or 



m co mpact groups of galaxie s (jlglesias-Paramo & Vflche? 



sjigle 

200JI; iMen des de Oliv eira et al j [2001; Temporin et al. 2003i: 
Lopez-Sanch ez et alj |2004; Am ram et al.l 2004; Lisenfe ld et al.l 
2004) while modelling has shown that the formation of conden- 
sations is possible in mergers b etween two disc g alaxies, as those 
happ ening in ULIRGs (see iDuc et al.l 120041 ; IWetzstein et al.l 
12005; Bournaud & Dudl2006l and others). 

In spite of the fact that they have the appropriate nature to 
harbour TDG candidates, no systematic search for these candi- 
dates among ULIRGs has been performed up to date . The TDG 
candidate in The Superantennae (Mirabel et al.l 19911) . turned out 
to be a background object and to our knowledge, only three more 
candidates have been identified so far (Mihos & Bothun 1998). 

An interesting question is whether these objects are likely 
to contain the same quantity or more TDGs than other types of 
merging systems. More generally, are ULIRGs favourable sites 
for star formation in and extended mode? Regarding the gas dis- 
tribution, it seems that ULIRGs are particularly efficient in driv- 
ing gas to the innermost regions, meaning that there is a priori 
much less gas reservoir in the outer parts. However, detailed sim- 
ulations able to explain the luminosities observed in these ob- 
jects and focused in the gas distribution show a dual behaviour 
for the gas: while the one in the inner disk (typically ~ 5 kpc) 
flows directly toward the central regions within 1 - 2 x 10^ yr 
after the pericent ric passage, the outer gas is ejected into tidal 
tails and bridges dlono et al.ll2004l). R egarding the extend of star 
formation in mergers, iBarnesI (120041) showed how this is larger 
when shock-induced star formation laws are used instead of lo- 
cal gas density based laws. From the observational point of view, 
the fact that ULIRGs are more gas-rich than less luminous in- 
teracting systems and the p resence of shocks in these systems 
dMonreal-Ibero et ani2006l) may favor extended star-formation. 
In this context, the positive or negative detection of TDGs (or 
external star formation in general) may help to clarify these is- 
sues. 

Finally, quantifying the incidence of TDGs in local ULIRGs 
may have cosmological implications since these objects are 
thought to be the local count erpart of the so-called Spitzer 
and sub-millimeter sources (e.g. Small et anil997l ; lHughes et alj 
1199 8) at z ~ 1-2. These sources p resent similar or even 
greater luminosities than ULIRGs (see iBlain et aril2002l for 
a review of their properties) and a relatively high fraction of 
them have morphological properties consistent with being sys- 
tems suffering an interaction or merging process (Chapman et alj 
2003; Popeetal. 2005; Small et al. 2004; Ivison et al. 2003 
Conselice et alj 120031) . They are two orders of magnitude 



more numerous th an local ULIRGs ( iPerez-Gonzalez et al.l2005t 
ICaputi et alj|2007l) . Thus, while the present-day ULIRGs could 
provide a relatively small contribution to the total number of 
newly created TDGs, their high-z counterpart might be funda- 
mental in this regard. Hence the detection of this kind of regions 
among ULIRGs may increase considerably the estimation of the 



percentage of today dwarfs that could have been formed from 
tidal debris. 

In this paper, we characterise the extranuclear star forming 
regions of a sample of ULIRGs using the combined information 
of Integral Field Spectroscopy (IFS) data together with high res- 
olution images from the HST. Properties such as extinction, ion- 
isation state, metallicity, age and mass of the stellar population, 
velocity dispersion, relative velocity etc., will be derived. We 
will use these parameters to estimate the likeliness of survival 
of these regions as future TDGs. The present work is part of a 
wider program whose final aim is performing a detailed study of 
a representative sample of ULIRGs using IFS. Due to the com- 
plex nature of these systems, this technique, that allows to obtain 
at the same time spectral and spatial information, is well suited 
for their characterisation. Previous re sults of this program can be 
found in lMonreal-Ibero et alj (|2006|) and references therein. 

Throughout the paper, a cosmology with 70 km s"' Mpc"', 
Qm = 0.3 and Q.\ - 0.7 is assumed. 

2. Sample and Observations 

2.1. Sample 

The galaxies under study were selected from a sa mple of nine 
ULIRG s with IFS data, including those presented in'Colin a et alj 
dlOOSl) plus IRAS 16007+3743 (Garci'a-Man'n 2005). They 
cover the ULIRGs low luminosity range (12.10 < log(LiR/L0) < 
12.60), have a variety of activity classes, and are in different 
stages of the interaction process. 

Potential TDG candidates have been selected as any high 
surface brightness compact region in the emission line maps (ob- 
tained from the IFS data) at a distance from the nucleus of the 
galaxy larger than 2'.'0, and associated with a single condensa- 
tion in the WFPC2 image. 

Following the above selection criteria, no regions are found 
in Mrk273 (Colina et al. 1999), Arp22q (Arrib as et al. 2001) and 
IRAS 17208-0014 jArribas & Colinal l2003h in the IFS data. 
IRAS 15206+3342 presents a chain of knots at a projected dis- 
tance of about 3.5 - 7.5 kpc f rom the nucleus in the WFPC2 
image (Arribas & Colina'2002'). However, the spatial resolution 
of the IFS data doesn't allow us to derive their individual prop- 
erties and, therefore, this system will not be considered here. 

Thus, the final sample include the following five sys- 
tems with at least one region of interest: IRAS 08572+3915, 
IRAS 121 12+0305, IRAS 14348-1447, IRAS 15250+3609 and 
IRAS 16007+3743. IRAS 15250+3609 is an advanced merger 
while the rest, with two nuclei and a more spread star formation, 
are in an earlier stage of the merging process. Table [1] shows 
some relevant properties of the selected galaxies. With the ex- 
ception of IRAS 16007+3743, which is at a slightly higher red- 
shift, typical redshifts for these systems are z ~ 0.07 which im- 
phes a linear scale of about 1.3 kpc arcsec 

2.2. Observations 

The IFS data we re ob tained with the INTEGRAL sys- 
tem JArri bas et al jl 19981) plus the WYFFOS spectrograph 
jBingham et al.ll 19941) in the 4.2 m WHT at the Observatorio del 
Roque de los Muchachos (Canary Islands) on April 1998 and 
April 2001. Spectra were taken using a 600 lines mm ' grating 
with an effective resolution of 4.8 A. We used the fibre bundle 
SB2 which has a field of view of 16'.'0 x 12'.'3 and is made up 
of 219 fibres, each 0'.'45 in radius. Fibres are arranged in two 
sets which observe simultaneously the target and the sky. The 
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Table 1. ULIRGs sample 



Galaxy 


z" 


Scale 




IR Class"^ 


Interaction 


Regions^ 






(kpc arcsec"') 






class."* 




IRAS 08572+3915 


0.058 


1.13 


12.17 


W 


m 


2 


IRAS 12112+0305 


0.073 


1.39 


12.37 


C 


m 


5 


IRAS 14348-1447 


0.083 


1.56 


12.40 


C 


m 


1 


IRAS 15250+3609 


0.055 


1.07 


12.09 


C 


IV 


1 


IRAS 16007+3743 


0.185 


3.10 


12.11 


C 


III 
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" Redshifts taken from the NASA/IPAC Extragalactic Database (NED). 

Infrared luminosities calculated using the infrared IRAS fluxes of lMoshir et alj ( I1993I) . the L/r expression given in I Sanders & Mirabell dl996t) 

and a luminosity distance assuming /fo = 70 km Mpc-\ CIm = 0.3 and Ha = 0.7. 
= Standard IR classification where warm (W) ULIRGs have /25//60 > 0.2 while cold (C) ULIRGs have /25//60 < 0.2. 

Following the criterion proposed by Veilleu x et al. (2002) where m means pre-merger and iv, merger state. 

Number of external star-forming regions analysed in this work. 



covered spectral rang e, exposure time and air ma ss for each ob- 
ject can be found in iMonreal-Ibero et aL I (12006 1). The data for 
IRAS 16007+3743 were obtained with a similar configuration 
and with a total exposure time of 6x1500 s. 



3. Data Reduction and analysis 

General reductions for the present IFS dat a have been discussed 
somewhere else (e.g Monr eal-Ibero et aTl l2006, and references 
therein). Here we will detail the procedure for the absolute cal- 
ibration of the data, which has not been included in previous 
papers in this series and is of particular relevance for deriving 
the Ha luminosity of the regions under study. 

3.1. INTEGRAL data absolute flux calibration 

IFS seems to be a suitable technique to carry out the absolute 
flux calibration in objects like ULIRGs, with a complex morpho- 
logical and kinematic al structure and at a certain redshift. Other 
options have greater difliculties. For example, observations with 
long-slit present centring problems, worsened in presence of dif- 
ferential atmospheric refraction. On the other hand, observations 
with narrow filters are complicated to calibrate due to redshift, 
internal movements of the system and blending of some lines, in 
particular Ha+ [N 11] . 

As a first step, the relative flux calibration was carried out 
using the fibre with greater S/N in the calibration star image 
to create a sensibility function which was used to calibrate ev- 
ery spectrum, correcting from possible effects of differential at- 
mospheric refraction ( Arribas et al. 1999; Filippenko 19 82), and 
from the light lost between fibres which affect the calibration 
star observations. Note that one can infer the fraction of the total 
flux collected by the selected fibre from images generated from 
the IFS data at different wavelengths (in a AA) . This informa- 
tion allows us to find the conversion factors counts/absolute flux 
and, therefore, calibrate absolutely the IFS data. In practise we 
proceed as follows. Firstly, we divided the calibration star spec- 
tra in several sections of a few Amstrongs each. Secondly, using 
the flux measured for each fibre, and its position in the bundle, 
we created a flux map of the calibration star for each section. 
Finally, we measured in that map the flux collected by the se- 
lected fibre and the total flux of the star in arbitrary units. The 
ratio between the flux collected by the fibre to the total flux of 
the star in the map allows us to convert the counts from the star 
into absolute flux in the selected wavelength. This is repeated 



for a set of images at different wavelengths, providing the con- 
version curve. In this specific case, we obtained the conversion 
factor generating images (at 0.04 arcsec pix"') in ten sections of 
300 A each through the entire observed spectral range. 

Three tests were performed in order to estimate the good- 
ness of this calibration method. Firstly, we compared the to- 
tal measured flux in a certain spectral range for an image of a 
star once calibrated with additional calibration stars available 
for the night. Secondly, the internal consistency of the calibra- 
tion method was checked by comparing measurements in Ha 
and Hj6 in certain emission line regions for one of our targets 
(IRAS 12112+0305) which had been observed in two differ- 
ent nights, pointings and rotation angles. Finally, the Ha flux 
measured in the nuclei of this system were compared with those 
found in the literature (Kim, Veilleux, & Sanders 1998) simulat- 
ing the aperture of an slit in the interpolated maps. From all these 
control tests we estimate uncertainties in the flux calibration of 
about 10 - 15%. 



3.2. HSTdata 



Data archival images from the HST taken with the WFPC2 
(F814W filter) were also used to complete our analysis. All im- 
ages but those of IRAS 08572+3915 were taken in snapshot 
mode and with the optional parameter CR- SPLIT activated. For 
every pointing there were two already reduced images which 
were combined to reject cosmic rays. Integrated counts for a 
given aperture were converted to ma gnitudes in the Vega system 
using the expressions given in iBaggett et al.l (l2002h . 

F814W filter is similar to the Johnson-Cousin I filter 
Transformations between the WFPC2 and the Johnson-Cousin 
filter sys tems need at least measur ements in two filters (see for 
example lOrigfia & Leithereij|200C ft. As in most of the regions, 
the only HST measurement available was that with this filter, we 
have decided not to apply any transformation and to use directly 
the F814W magnitudes as if they were I magnitudes. We have 
performed a comparison between the expected magnitudes for 
the F814W and I filter with SYNPHOT using different synthetic 
spectra. The more relevant for this study was the one correspond- 
ing to an instantaneous burst of 6 Myr with different levels of 
extinction. Differences were always 0.1 mag or less. 
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4. Results 

4.1. Identification of the regions of interest 

The selected regions are identified as those bright condensa- 
tions in the emission Unes outside the nuclear region (i.e. at dis- 
tances of at least 2 kpc from the nucleus). They are marked in 
the first column of Figure [T] which contains the WFPC2/HST 
images. The other two columns in this figure show the contin- 
uum and Ha emission line maps. The stellar continuum was 
obtained averaging two continuum ranges towards the blue and 
red of the Ha emission line simulating the action of a filter. For 
the ionised gas component, we have fitted Gaussian functions 
to the emission lines using DIPSO package (iHowarth & Murravl 
119881) within the STARLINK environment (details can be found 



Monreal-Ibero et 31112006 '). The general morphology of the 



m 

stellar and ionised gas component are similar although the stel- 
lar component presents a more spread distribution than that of 
the ionised gas. In the following, we describe briefly the identi- 
fied regions: 

IRAS 08572+3915: Two condensations have been identified 
here: k7 at ~ 6" from t he northern g alaxy toward the north- 
east, neai- its ti dal tail ([A rribas et alJuOOQ) and Rl (k4+5 in 
the notation of I Arribas et al. 2000^in the southern tidal tail 
and at a distance of ~ 5'.'7 from the southern nucleus. Here 
k4+5 will be treated as a single region as we were not able to 
close independent contours for them in the WFPC2 image. 

IRAS 12112+0305: This system shows several regions of inter- 
est. The brightest emission line region (Rl) is not associated 
with any of the nuclei of the system, but with a low surface 
brightness r egion in the WFPC 2 image towards the east of 
the system dCoUna et al.l l200dl) . This image also shows an 
arc-like condensation located 3" north of the southern nu- 
cleus at PA. ~ 20° (R2) which is relatively bright in the Ha 
emission line map. Finally, the southern tidal tail shows a se- 
ries of condensations which correspond either to peaks in the 
emission line map (kl and k2) or in the continuum map (kc). 

IRAS 14348-1447: The most interesting region in this system 
has already been identified by Mihos & Bothun ( 1998) as a 
peak of Ha emission towards the southeast of the system (R). 
In addition, there are other two regions which should be men- 
tioned here. On the one hand, the peak of Ha emission in the 
northern galaxy is not associated with its nucleus but with 
four knots located at the base of its tidal tail. Similarly to 
IRAS 15206-1-3342, the individual knots cannot be resolved 
in the ground based data and its analysis will not be con- 
sidered here. On the other hand, the WFPC2 image shows a 
faint condensation at ~4'.'0 and P.A.~ 70° from the southern 
nucleus whose analysis will also not be considered here due 
to the elevated flux contamination from the southern galaxy 
in the INTEGRAL data. 

IRAS 15250+3609: This galaxy presents a low surface bright- 
ness region in the F814W band at ~7'.'0 from the nucleus to- 
wards the southwest (R), w hich is however relativ ely bright 
in the line emission maps dMonreal-Iberol l2004h . HST in- 
frared images show a series of knots in the main body 
of the galaxy probably associated with globular clusters 
dScoville et al]|2000 ). Although they can be responsible for 
some of the extended emission observed in the ground based 
data, they cannot be identified as individual sources and their 
analysis is not considered here. 

IRAS 16007+3743: The morphology of IRAS 16007+3743 
could indicate that this sy stem is a multiple merger dCui et alj 
l200UlBorne et alj|2000l) . However, kinematical information 



derived from INTEGRAL (i.e. velocity field and velocity 
dispersion maps) made us classify this system as a system 
of two spiral galaxies in interaction ( Garci'a-Marin 2005'). 
Three regions have been identified: Rl and R2 seen in the 
Ha map as two bright condensation along the northern tidal 
tail and R3, a region at ~3'.'0 towards the south of the system 
apparently not associated with any of the tidal tails. 

To have an idea of the data quality, the observed spectra for 
the regions under study are shown in Figure |2l From blue to 
red, Hfi, [O iii]^^4959,5007, [O i]^6300, Ha+[N ii]^^6548,6584 
and [S ii]/l/l67 17,6731 are detected in most of the regions with 
enough signal-to-noise. 

4.2. Cliaracterisation of tlie external star-forming regions 
4.2.1 . Line ratios 

Measured line ratios are shown in Table |2] and Figure 
3] which contain the cla ssical diagnostic diagrams of 



Veifleux & Osterbrockl (11987 ). The regions present line ratios 
typical of H n regions. Only by means of the [O i]A63QQ/Ha line 
ratio some of them could marginally be classified as LINER. 
Star formation is the most plausible mechanism to explain the 
observed line ratios i.e. the line ratio of all the regions analysed 
here are well within the limits proposed bv lKewlev et alj d2001h 
for ionisation by star formation. 

In general, lin e ratios are comparable to those found in TDGs 
can didates (e.g. | Duc & Mirabe I ll998t iTemporin et al] 120031: 
Lop ez-Sanchez et al.1 |2004|) . except for flieir [Oi]A6300fHa 
which is smaller by ~0.3 dex in those cases where it could be 
measured. However, the poor statistics (i.e. the [Oi]/I6300 line 
could only be measured for four regions) makes difficult to as- 
sess the significance of this result. 



4.2.2. Metallicities 

One observational characteristic that helps to establish the 
tidal origin of a certain dwarf galaxy is its high metallic- 
ity, whic h suggests that TDGs are made up from processed 
material dPuc & Mirabe]|[T998l) . Abundances are usually esti- 
mated using empirical methods based on the intensities of cer- 
tain optical lines like the w idely used m ethod of the R23 cal- 
ibrator (i.e. [j orres-Peimber t et alj[r989l) and the 5 23 indicator 
dVflchez & Estebaal996) . However, both of them involve emis- 
sion lines which are not available within our spectral range. 
Nevertheless, it is possible to estimat e the m etallicity using the 
N2 calibrator proposed by Denicolo et alj d2002) - based on 
the ratio between the [Nn1/ I6584 and Ha emi s sion li nes -, and 
the empirical diagrams of lEdmunds & Pagel dl984l) . - based 
on the [ O iii]/l/l4959,5007 to H /? ratio - parametrised as ex- 
plained in lDuc & Mirabell dl998'). The predicted metallicity val- 
ues expected from these indicators are shown in the last two 
columns of Table |2] Both indicators predict similar values for 
the metalHcity (average of 12 + logO/H - 8.70 for the N2 
calibrator against 8.57 for the one involving the oxygen lines). 
Giving that typical uncertainties for the line ratios are about 
0.2 dex (which translates into ~0.3 for the metallicity), differ- 
ences between both indicators are not significative. This rela- 
tively high metallicities are about a factor 2 larger than those de- 
rived for other TDGs candidates (-8.35, Due & Mirabel 199|; 
IWeilbacher. Due. & Fritze-v. Alvensleben[i2003h and about half 
the solar metallicity. 
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Region 






E(B - Vf 




[Niil/Ho''''^^ 


[S iij/Ha^" 


12 + log(0/H)'' 


12 + log(0/H)'= 




(") 


(") 










([OlIl]/^8„ppe,-) 


(N2) 




IRAS 08572+3915* 


Ris 


5.2 


-7.4 




0.24 


-0.85 




8.49 


8.50 


k7 


5.3 


2.2 




0.54 


-0.82 




8.29 


8.52 


IRAS 12112+0305 


Rl 


2.7 


-1.5 


0.25 


0.22 -1.24 


-0.66 


-0.58 


8.50 


8.64 


R2 


-1.5 


0.6 


1.07 


0.01 -1.02 


-0.41 


-0.43 


8.65 


8.82 


kc 


-1.8 


-7.4 


0.47 


0.16 


-0.64 


-0.99 


8.55 


8.65 


kl 


-2.8 


-7.2 


0.66 


-0.12 


-0.49 


-0.38 


8.74 


8.76 


k2 


0.7 


-7.7 


0.45 


0.08 


-0.42 


-0.39 


8.60 


8.81 


IRAS 14348-1447 


R 


5.1 


-1.7 


1.12 


0.26 -1.36 


-0.55 


-0.49 


8.48 


8.72 


IRAS 15250+3609 


R 


-5.0 


-4.8 


0.25 


0.26 


-0.78 


-0.60 


8.48 


8.55 


IRAS 16007+3743 


Rl 


0.9 


5.4 


0.55 


0.13 -1.05 


-0.47 


-0.67 


8.59 


8.78 


R2 


0.7 


2.8 


0.40 


-0.02 


-0.55 


-0.67 


8.65 


8.71 


R3 


2.2 


-1.9 


0.75 


0.49 


-0.55 




8.59 


8.75 



^ Relative positions of the regions. We took the northern nucleus as reference in IRAS 08572+3915 and IRAS 12112+0305; the southern one 
in IRAS 14348-1447 and the western one in IRAS 16007+3743. 

Calculated as the average value in a O'.'45-radius aperture centred in the region of interest. 

Extinction corrected data. We assumed a foreground screen model, and used the reddening curve of lWhitfor? ('1958) parametrised as explained 

in Miller & Mathews ( 1972) and R=3.1 (Rieke & Lebofsky 1985), except for regions Rl and k7 of IRAS 08572+3915 where no extinction 

measurements are available. Typical errors for line ratios are <0.2 dex. 

Derived metallicities using the upper branch of the empirical d i agram of lEdmunds & Pagel ( 119841) . 
° Derived metallicities using t he N2 calibrator of | penicol6 et alj ( l2002l) . 
' Measured values taken from lArribas et alj ( |2000|) . ^ 

^ This corresponds to knots k4 and k5 in the notation of lArribas et alj ( I2OOOI) . Their relative positions are: (5.8,-7.8) and (4.7,-7.0). 



4.2.3. Sizes 



4.2.4. Ha equivalent wicJtlis and luminosities 



Sizes were estimated from the WFPC2/F814W images. The ir- 
regular shape of the regions as well as some possible contamina- 
tion from other structures within the system (i.e. other regions, 
tails or the parent galaxy itself) made difficult to define their lim- 
its. In order to proceed in a systematic way and compare with 
the sizes of other extragalactic objects, we defined the size of 
a given region as the area encircled in the largest closed con- 
tour centred in the emission peak associated with that region. 
This allows us to define an equivalent radius as rgqu - -^Aiiea/n 
which gives an estimate of the total size of the region. An effec- 
tive radius was also measured as the one which contains half of 
the flux within this area (see columns 7 and 8 of table [3]). They 
can be use for comparison with the characteristic radii derived 
for other objects. The ratio between the effective and equivalent 
radii (0.3 - 0.7) gives and idea of the compactness of a certain 
condensation. 

Equivalent radii range from a few hundreds pes to ~ 1.5 kpc 
(mean of 750 pc). In general, when comparing with the sizes 
of Hii regions (which have typical radii of ~ 100 - 900 pc), 
the p resent regions are similar to the largest Giant H 11 regions 
(lKennicutt.1984; .Mavva..l9 94). 



All the regions under study have effective radii (mean ~ 
430 pc) comparable to dwarf galax ies located in the Local 
Group (r ~ 0.3 kpc, iMateol (l998 f). to the so-called 

Cnmnart Dwarf tralaxips with effect\vp radii ?<r„«<l 



Blue 



(e.s.lMarlowe et al. 


1997; Cairos et al. 2003) or to some already 


detected TDG (e.g. 


Duc&Mu-abel,1998,). 



Ha equivalent width and luminosity distributions are shown in 
Figure |4] together with the values measured for some TDGs 
as well as Hii regions taken from the literature. The rela- 
tively high Ha equivalent widths derived (i.e. 34 - 257 A), 
indicative of the existence of a very young stellar population, 
are comparable to those measured for TDG candidates (e.g. 
12 - 623 A forllglesias-Paramo & Vflchezll200TI or 24 - 198 A 
for 'Temporin et al.' '2003") and for extragalactic Hn regions in 
general (Mavva 1994). Regarding the Ha luminosities, once 
the extinction effects have been taken into account, the se- 
lected regions are more lumino us than the Giant H ii Regions 
(lMavvalll994l; lKennicuttlll984l) . In general, these luminosities 
are in better agreement with those found for TDG candidates 
(Due & Mirabel 1998; Iglesias-Paramo & Vilchez 2001). The 
most suitable sam ple for com parison with our candidates is the 
one of 1 Temporin et al.l (l2003i) which takes into account extinc- 
tion effects. With the exception of three very bright regions (see 
discussion below), the Ha luminosities are comparable to those 
of Temporin et al. (mean ~ 3 x lO**" erg s"' for their sample 
against ~ 6 x 10"*" erg s"' for ours). 

The extinction corrected luminosity of three regions are well 
above 5 x 10"*' erg s"' which is more than one order of mag- 
nitude larger than the luminosities found in typical TDGs. Two 
of them belong to IRAS 16007+3743 which is at z=0.185 and 
has a factor of 2 - 3 lower linear resolution than the rest of the 
sample. One possibility is that emission from these two particu- 
lar regions was due to a wider area that encompasses more than 
one condensation. However, a resolution effect cannot be used 
as explanation for the high luminosity of the third region (R2 
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Table 3. Properties of the Star-Forming Regions (II): Ha related observables, I magnitudes, sizes and distances. 



Region 


EW(Hq')'' 
(A) 


(10"'* erg cm-2) 


L(Haf itif 
(10^° ergs-') 




Tequ 

(pc) 


fcff 

(pc) 


(kpc) 


(kpc) 




IRAS 08572+3915" 


Rl 

k7 


78 
245 


7.0 
4.0 


0.6 23.78 
0.3 22.62 


-13.30 
-14.46 


373 
239 


121 
128 


8.8 
6.4 


4.2 
6.4 


IRAS 12112+0305 


Rl 
R2 
kc 
kl 

k2 


257 
79 
34 
46 
71 


31.8 
37.4 
2.7 
3.5 
3.4 


7.4 21.12 
57.6 19.05 

1.0 21.15 

2.1 21.43 

1.2 22.02 


-16.86 
-20.14 
-17.16 
-17.16 
-16.26 


373 
698 
451 
383 
378 


224 
494 
311 
150 
249 


5.0 
2.4 
9.0 
9.0 
9.0 


4.3 
2.2 
6.9 
7.1 
8.4 


IRAS 14348-1447 


R 


194 


15.5 


34.1 20.40 


-19.13 


1082 


546 


8.2 


6.0 


IRAS 15250+3609 


R 


141 


6.5 


0.6 20.73 


-16.36 


891 


420 


6.9 


6.9 


IRAS 16007+3743 


Rl 
R2 
R3 


40 
234 
80 


12.0 
35.5 
5.85 


57.1 20.16 
189.5 19.80 
8.8 21.83 


-20.63 
-21.07 
-19.23 


1151 

1375 
1562 


828 
884 
851 


17.1 
9.4 
5.8 


16.9 
8.9 
9.0 



^ Calculated as the average value in an O'.'45-radius aperture centred in the region of interest. 
As a compromise between collecting all the flux of the region and avoiding the contamination from the neighbouring zones of the systems, 
flux was measured in a I'.'O-radius aperture in all regions of IRAS 08572+3915, IRAS 14348-1447, IRAS 15250+3609 and Rl and R2 of 
IRAS 121 12+0305 while for those in the tidal tail of this system and those in IRAS 16007+3743, a O'.'5-radius aperture was used. 

° Extinction corrected Ha luminosity except for regions Rl and k7 of IRAS 08572+3915 where no extinction measurements are available. 
Observed I magnitude, measured within the area enclosed by the biggest closed contour centred in the emission peak (see explanation in text). 

^ Extinction corrected absolute I magnitude. 

' Distance to the mass centre of the system. 

^ Distance to the nearest galaxy. 

'' Ha equivalent widths and fluxes taken from lArribas et alj d2000h . 



of IRAS 12112+0305, scale -1.4 kpc arcsec"'). Alternatively, it 
may be possible that the relatively large amount of gas involved 
in ULIRGs in general, and the violent event suggested by the 
morphology of IRAS 16007+3743 and IRAS 121 12+0305 could 
lead to the large luminous TDGs. 



It is widely known that Ha emission can be used to con- 
strain the properties of recent episodes of star formation. An 
upper limit to the age can be estimated by using its equivalent 
width, as it decreases with tim e (see inset in Figure [51 ). We have 
used STARBURST99 models dLeitherer et al.|[l999h to estimate 
ages for the bursts in the extranuclear regions. In view of the 
estimated metallicities in section 14.2.21 we chose two different 
spectral synthesis models, both for an instantaneous burst with 
a Salpeter IMF, and upper mass limit of 125 M© but with two 
different metallicities: Z = 0.020 and Z = 0.008. Models with 
continuous star formation cannot explain the observed equiva- 
lent widths as the newly born massive stars are able to maintain 
too elevated values of the equivalent width in comparison to the 
observed ones even for ages greater than 100 Myrs. The esti- 
mated ages for the bursts are shown in Figure |5] and second col- 
umn of tablelD They range between 5-8 Myr, depending on the 
assumed model, confirming the presence of a young population. 



Summarising, equivalent widths of the analysed regions are 
similar to those derived for both extragalactic H n regions and 
TDG candidates and they are typical of young burst of star for- 
mation. The Ha luminosity place these regions more in the group 
of TDG candidates than in that of extragalactic H ii regions. 



4.2.5. I magnitudes 

We have measured the I magnitudes within the apertures defined 
in section |4!23] using the task POLYPHOT under the IRAF envi- 
ronment. Typical errors in the measurements where ~ 0.1, with 
the exception of Rl in IRAS 08572+3915 that was ~ 0.4. The 
comparison between these magnitudes and those derived from 
STARBURST99 are shown in Figure |6] where plotted magni- 
tudes were obtained as an average of the predicted magnitudes 
for the Z = 0.020 and Z - 0.008 instantaneous burst models. 
The trend is good although the predicted values from extinction 
corrected Ha luminosities overestimate the observed magnitudes 
by ~0.9 mag on average. These differences can be explained by 
both, the uncertainties associated with the selection of the aper- 
ture in the different bands and the existence of other stellar pop- 
ulations either from a previous burst during the merging process 
or from the original population in the merging galaxies, as found 
in some TDG candidates (e.g. lLopez-Sanchez et alj|2004 ). The 
presence of such a older population would cause an underesti- 
mation of the equivalent width and thus, overestimation of the 
age of the burst. Younger ages for the bursts translates into a de- 
crease of the predicted I magnitudes of ~ 0.5-1.0 mag according 
to the models. 

4.2.6. Distances 

The two last columns in Table[3]show the projected distances be- 
tween the studied regions and the mass centre of the system, as 
well as the closer galaxy. Typical measured distances are ~8 kpc 
with a maximum of ~17 kpc. However, these objects were not 
selected based on a particular favourable geometry. Indeed, their 



A. Monreal-Ibero et al.: TDG candidates in ULIRGs 7 
Table 4. Properties of the Star-Forming Regions (III): Derived characteristics of the stellar populations and dynamical parameters. 



Region Age" m^^ a" M^y, M^/ M^J 

(lO^Myr) (10'^ Mo) (kms-') (10** Mq) (10** Mq) (10** Mq) (kms"') (kms"') 



IRAS 08572+3915 



Rl 


5.4 


1.9 




0.04 


0.09 


196 


133 




k7 


6.4 


2.2 


47 


5.9 0.03 


0.02 


41 


156 


196 


IRAS 12112+0305 


Rl 


5.3 


22.8 


60 


16.8 1.22 


1.10 


75 


410 


14 


R2 


6.4 


420.6 


56 


32.4 72.32 


53.65 


-280 


592 


0.4 


kc 


7.3 


11.5 




0.37 


0.32 


-338 


306 




kl 


6.9 


19.6 




0.23 


0.21 


-207 


306 




k2 


6.6 


9.4 




0.18 


0.11 


-343 


306 




IRAS 14348-1447 


R 


5.7 


136.6 


50 


28.5 10.6 


14.78 


144 


402 


1.9 


IRAS 15250+3609 


R 


6.5 


2.8 


50 


21.9 3.62 


3.62 


170 


264 


6.0 


IRAS 16007+3743 


Rl 


7.1 


581.6 


61 


64.4 1.7 


1.1 


-41 


302 


38 


R2 


5.4 


678.3 


76 


106.8 17.7 


112.39 


92 


408 


0.9 


R3 


6.4 


62.8 


85 


128.5 106.65 


12.58 


338 


519 


1.2 



" The age and mass in young stars obtained as the average from the predictions for the low and high metallicity models. 
Calculated as the average value in an O'.'45-radius aperture centred in the region of interest. 
Tidal mass assuming the potential is created by a point mass in the mass centre of the system. 
Tidal mass assuming the potential is created by the nearest galaxy. 
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Fig. 4. Ha luminosities vs. equivalent widths. Data from the 
present study follow the same colour/symbol code as in Figure 
[3] We also include for co mparison mea surements for the ex- 
tragalactic Hii regions of iMavval d 19941) (dots, not extinction 
corrected) as well as for other extranuclear star forming re- 
gions which have been ana lysed as possible TDG candidates 
(solid symbol s). Those of Iglesias-Paramo & Vflche3 (1200 lb 
(squares) and iMundell et al. (12004') (hex agon) a re not extinc- 
tion corrected while those of Tempori n et alj ( l2003h (trian- 
gles) have been corrected for extinction. The vertical dashed 
line shows the luminosity limit for TDG candidates used in 
llglesias-Paramo & Vflche"3(l200ll) . 
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Fig. 5. Measured Ha equivalent widths in comparison to two 
SB99 instantaneous bursts of Z = 0.020 (continuum line) and 
Z = 0.008 (dashed line). Colour/symbol code is the same as in 
figure |3] The predicted ages range from 5.5 to 7.1 Myr if the 
model with solar metallicity is assumed or from 5.1 to 7.8 Myr 
if the one for lower metallicity is used. The small inset shows 
the covered range in equivalent width for both continuous and 
instantaneous STARBURSTS99 models (thin and thick lines re- 
spectively). Note the elevated equivalent widths predicted for the 
continuous star formation models, which are not reached by any 
of the regions analysed here. 



geometries are very complex and the influence of projection 
effects - which should be present - is particularly difficult to 



quantify. These distances have to be seen as lower limits to the 
actual ones. These are relatively close distances to the parent 
galaxy when compared with other condensations in Ha identi- 
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Fig. 6. Comparison between measured and predicted I magni- 
tudes. Colour/symbol code is the same as in Figure|3] The dotted 
Hue represents slope one. Predicted magnitudes were obtained 
as an average of the predicted magnitudes for the Z = 0.020 and 
Z - 0.008 instantaneous burst models. 



fied as TDGs (e.g. iDuc & Mirabelll 19981) although there are al- 
ready some examples of TDG candidates at such distanc es (see 
for instance TD4 4a in llglesias-Paramo & Vilchezll200ll or ob- 
jects A2 and E in lAmram et al ■1 12004 . 



5. Discussion: What are the chances of survival of 
the analysed regions as TDGs? 

Most of the observational properties of the selected regions (i.e. 
Ha equivalent widths and luminosities, metallicities, ionisation 
states and radii) derived in section |4] place these regions in the 
locus of the most luminous extragalactic Giant H ii Regions and 
are consistent with those expected for TDGs or TDG progeni- 
tors. In the next section, we will use the kinematical information 
provided by INTEGRAL which, together with the evolutive state 
of the system, will allow us to estimate the likeliness of survival 
of these candidates. 

As TDGs are stable entities with their own dynamic, the 
best defi nition for this kind of objects is the one proposed by 
IDuc et a l. (2000) who identify a TDG as that object that con- 
stitutes a self-gravitating entity and that was made up from the 
debris of a galaxy interaction. To evaluate the chance of survival 
of TDG candidates, two basic questions need to be answered: i) 
Is the candidate massive enough to survive to its internal move- 
ments? ii) Is the candidate massive enough to survive to the grav- 
itational forces exerted by the parent galaxy? Depending on the 
available observables, several criteria have been used in the liter- 
ature in the past to classify a certain condensation as a TDG. In 
the following we apply them (when possible) to our star-forming 
regions. 

5.1. Stability against internai motions 

esias- Paramo & Vflched ( l200lh established a luminosity cri- 
terion (L(Hq') > 10^^ erg s"') which should fulfil systems sta- 
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Fig. 7. Velocity dispersion vs. estimated eff'ective radius. The 
colour/symbol code is the same as in figure [3] The typical size 
for errors is shown in the lower right corner. The continuous 
line represents the fit for ex tragalactic H ii regions obtained by 
iTerlevich & Melnickl (Il98ll) while the dashed and dotted lines 
show their fit for elliptical galaxies, and globular cluster + el- 
liptical galaxies respectively. Small symbols represent different 
samples of d ynamically hot sys t ems: open circles, ga lactic glob- 
ular clusters jTrager et al.ll993l:lPrvor & Me vlan'199 3i); crosses , 
massive globular clusters in NGC 5128"lMartini & Hol[2004 : 
open triangles, dwarf elliptical galaxies (iGeha et al.ll2003l) : open 
squares, intermediate ellipticals; solid circles, giant ellipticals; 
solid triangles, compact elliptical s; asterisks, dwarf ellipticals; 
and stars, bulges (all of them from lBender et alJll992l) . 



ble to internal motions. As it has been mentioned above (section 
|4.2.4| i all the selected candidates meet well this criterion (see 
also Figure 4), especially two very luminous regions of IRAS 
16007 -f 3743 and R2 in IRAS 12112-1-0305. 

One of the best test to check if a certain detected con- 
densation constitutes a self-gravitating entity is detecting ve- 
locity gradients that could indicate the e xistence of indepen- 
dent rotation (e.g.LW eilbacher et al. 2002; Bournaud et alj|2004l; 
iMendes de Oliveira et al..2001,) . However our data lacks the nec- 
essary angular resolution to resolve any velocity field across the 
extranuclear condensations. Typical size of the detected regions 
is roughly equivalent to 1-3 SB2 fibres, which doesn't allow to 
derive the individual velocity field of the regions under study. 

Another method that can be used to establish whether or 
not the TDG candidates are stable, self-gravitating entities is 
to study their location in the radius-velocity dispersion cor- 
relation measured for ellipticals and globular clusters. With 
this purpose, we have plotted the classical fits in the r^q - cr 
plane for elliptical galaxies, globular clusters and Hii regions 
of Terlevich & Melnick ( 1981) together with our data in Figure 
|2] For the present INTEGRAL configuration the spectral res- 
olution is cr ~80 km s We considered that we were able 
to determine a reliable measurement of the velocity dispersion 
when this was greater that half the spectral resolution. We have 
also overplotted the position of recent samples of hot dynam- 
ically system such as galact ic globular clusters dTrager et al.l 
11991 lPrvor(feMevlanl[l993h . massive globular clusters in in 
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NGC5128 dMartini & H o 2004), dwarf elliptical galaxies in the 
Virgo cluster 'Geha et al. ( 200 3p. ari d the sample of dynamically 
hot systems of Bender et al] ( |1992|) . Based on their position in 
the radius-velocity dispersion diagram, the selected TDG candi- 
dates should be self-gravitating entities. Moreover, they are in a 
region close to other dwarf gala xies (compare in particular with 
the sample of iGeha et aLll2003l) . and far from the region occu- 
pied by globular clusters. 

5.2. Stability against forces from tine parent gaiaxy 

One way to see whether a condensation is stable against the grav- 
itational p otential of the parent galaxy/system is estimating its 
tidal mass (iBinnev & Tremain^ll987l;lMendes de Oliveira et alj 
1200 lb which is defined as: 

M„,, = 3m(^) (1) 

where M is the mass of the parent galaxy, R is the radius of 
the TDG candidate (here estimated as r^g,,) and D is the distance 
to the parent galaxy. If the tidal mass is smaller than the mass 
for the TDG candidate, then the candidate is stable against the 
forces exercised by the parent galaxy. This expression is valid in 
those cases where the size of a certain regio n is small in compar- 
ison w ith the distance to the parent galaxy jBinnev & Tremain^ 
Il987 l). Typical values for the ratio size / projected distance 
(which is an upper limit for the real size / distance) for our con- 
densations range between 0.04 (k2 in IRAS 12112-1-0305) and 
0.32 (R2 in IRAS 12112+0305), small enough to try to asses 
their stability against forces from the parent galaxy by means of 
the use of the tidal mass. 

In general, the gravitational potential for a ULIRG is a 
complex function of the mass distribution of the system which 
evolves with time but, as a first approximation, it can be consid- 
ered dominated by the masses of the main bodies of the inter- 
acting galaxies. We will consider here two simple approaches: 
assuming that the TDG candidate is suffering the gravitational 
potential i) of a point mass in the mass centre and with the total 
mass of the system (M^^) and ii) of the nearest galaxy (M^j^). 
The estimated values are shown in table |4] For Mf.^ we used 
M — total mass of the system and D — distance to t he mass cen- 
tre. M asses for the parent galaxies were taken from lCohna et al.l 
( I2005h and iGarcfa-Marfnl (l2005h . For the final tidal mass (M,,^) 
we have considered the most conservative approach: the maxi- 
mum of M'^.j and M^.^. It has to be taken into account that we are 
measuring the projected distance (always smaller than the real 
distance). Therefore, these values have to be seen as an upper 
limit to the real tidal mass. 

An estimation of the dynamical mass of the candidates is 
needed to compare with the tidal masses. For those regions 
where a reliable measurement of the velocity dispersion was 
achieved, this can be done using: 

M{Mq) = cte X 10''/?eff(kpc)cr2(km s"') (2) 

Estimations of the constant range between 1 .4 dBender et alj 

11992) and 2.2 (Terlevich & Melnick 198ll) depending on the 
method used for its determination and the mass distribution 
assumed. Here have adopted cte - 2.09, a value derived 
from the Virial theorem, assuming a de Vaucouleurs profile 
(which describes the luminosity distribution of spheroidal sys- 
tems quJte_well}jthat^ traces the mass of the sys- 
tem dCombes et alj [20021) and that the tridimensional velocity 
dispersion is related with the observed velocity dispersion as 



(y2) = 30-2 dMcCradv et al.ll2003h . As it can be seen in last col- 
umn of table 4, most of the candidates with an estimation of the 
dynamical mass have a MdxnlM,id ratio which supports the idea 
that most of the regions are stable against the forces exerted by 
the parent galaxy (mean -30). Only Rl in IRAS 12112+0305 
and R2 in IRAS 16007+3743 do not fulfil the criterion, but tak- 
ing into account the uncertainties involved this is not very sig- 
nificant even for these two cases. 

A last test could be a comparison between the relative ve- 
locity of a certain region and its escape velocity. This criterion 
cannot be blindly used. Firstly, it is a function of the distance 
which will be always larger than the projected one. Secondly, 
only one component of the velocity is being measured without 
having any information about the movements on the plane of the 
sky. Finally, as for a certain configuration it is not possible to de- 
termine if the region is closer or further from the observer than 
the mass centre, two possibilities for the relative movements be- 
tween the region and the system are always possible. Due to all 
these uncertainties we cannot give to this comparison the same 
importance (weight) as the comparison between the dynamical 
and tidal masses. However, we have found useful to include it 
here for completeness as used together with other indicators can 
give us a better idea about the survival chances of the candidates. 

In a simplified way, we consider that the gravitational poten- 
tial is created by a point mass in the mass centre and with the 
total mass of the system. Relative and escape velocities can be 
compared looking at columns 7 and 8 of table |4l As average, the 
absolute values of the relative velocities are much smaller than 
predicted escape velocities (189 km s"' vs 346 km s"^). In or- 
der to take into account projection effects, we will consider that a 
certain region will have good survival chances following this test 
if it satisfy the condition \vrei\ - |vv.se x cos(;7r/4)| > 0. According 
to this comparison, three out of twelve regions would pass this 
last test. 

5.3. Comparison with modeis 

Although the observational characteristics of these condensa- 
tions give some hints that make us to consider them as puta- 
tive TDG progenitors, a comparison with dynamical evolution- 
ary models is needed in order to better predict if they will be- 
come independent objects. 

The formation of TDG inside the tidal tails of inter- 
acting galaxies has b een the subject of recent simulations. 
IWetzstein et al.l ( l2005h showed the crucial role of the gas in 
the formation of these systems, contrary to the findings of 
iBarnes & HernquistI (ll992E . lBournaud & Du3 (12006) focused in 
both their formation and evolution of this kind of condensations. 
Typical masses of the surviving ones are greater than 10*^ Mq. 
These simulations end at ~2 Gyr after the beginning of the 
encounter As an average ~ 3.2 condensations per galaxy are 
formed. From those, only the most massive (>10^ Mq) and born 
at the tip of the tidal tail (~ 0.6 - 1.0 Rtaii) are found to survive 
to the encounter 2 Gyr after the pericenter We have performed a 
comparison of the properties derived for our condensations with 
those created in these simulations. Seven out of our twelve re- 
gions have more than lO'' Mq. Attending to Bournaud & Du3 
(120061) . these have more chance to survive to the interaction pro- 
cess, loosing only a small fraction of their mass. If the fate of 
our regions were similar to the one of these condensations, most 
of them (with the exception of R in IRAS 15250+0305 and Rl 
in IRAS 16007+3743) would probably fall back to the parent 
galaxy loosing part of its mass. A summary of the comparison 
with these simulations is shown in the sixth column of table |5] 
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One question that arise after this comparison is what is the ori- 
gin of regions such as R in IRAS 14348-1447. Its observational 
properties indicate that it is a condensation of matter made out 
of the material of the interaction which is suffering a burst of 
star formation and hence, it should, in principle receive the same 
treatment as the other regions as possible entities able to sur- 
vive to the merging process. However, to our knowledge, there 
is no simulation able to predict condensations of matter, even 
transient, outside the tidal tails. 

How will their evolution be after 2 Gyrs? Unfortunately 
there are not evolutionary models coverin g this ph ase, for the 
range of masses of ou r regions. Howeven lKroupal (1 19971) and 
iMetz&Kroiipal ( 120071) simulated the evolution of condensations 
with an initial mass of 10^ Mq (i.e. more than one order of mag- 
nitude lower than the present case) up to 10 Gyrs. They show that 
around 4-6 Gyr, depending on the model, the condensations re- 
tain ~ 10 % of their initial mass, and only ~ 1 % at 8-10 Gyr. If 
these results could be scaled to our more massive regions, they 
would end up with masses of about 10^ - 10^ Mq. However, this 
extrapolation may be questionable and a more firm conclusion 
should wait for reahstic simulations for the appropriate range of 
mass and time. 

5.4. Summary 

Due either to projection effects or to observational constrains, 
none of the tests performed above can alone tell us whether a cer- 
tain region is going to survive as a future TDG or not. However, a 
region that satisfy most of the tests has more chances to survive. 

We have summarised the results of these tests in table |5] We 
have assigned a High probability of survival to those regions that 
satisfy all the utiHsed criteria (i.e. only R of IRAS 15250+3609 
andRl in IRAS 16001 +31 43), High-Medium and Medium pmh- 
ability to those that satisfy four and three of the utilised crite- 
ria respectively, and Low probability to those that satisfied only 
two or one. The most promising candidate is the one detected 
in IRAS 15250+3609. The relative large projected distance to 
the parent galaxy (6.9 kpc), the advanced state of the merging 
process, its position in the r - cr plane and the large velocity 
difference with relati on to the parent galaxy dColina et alj|2005t 
lMonreal-Iberol2004b make it the best candidate to survive to the 
merging process. This region shows the way to be followed in 
the hunting of such systems: a systematic search for regions with 
recent or ongoing star-formation among a sample of ULIRGs 
made out of galaxies in an advanced state of the merging pro- 
cess seems appealing. The other region that should be mentioned 
here is Rl of IRAS 16007+3743. Appai-ently it is at the tip of 
the tidal tail of a system that is starting its interaction and de- 
pending of its evolution it could finish in a region similar to the 
one in IRAS 15250+3609. There are in addition three candi- 
dates with High-Medium probabiUty (Rl of IRAS 12112+0305, 
R of IRAS 14348-1447 and R3 in IRAS 16007+3743). Two of 
them are new candidates while R in IR AS 1448-1447 had al - 
ready been proposed as candidate bv IMihos & BothunI (Il998l) . 
The continuity in the velocity field between them and their par- 
ent galaxy reduce the chance of survival. 

The analysed systems are at higher redshift (and hence worst 
linear resolution) than their less luminous relatives previously 
studied. A follow-up of these candidates with IFS systems able 
to provide greater spatial resolution (i.e. Adaptive Optics as- 
sisted systems) would permit to look for in dependent veloc- 
ity gra di ents in them similar to the works of IWeilbacher et all 
(120021) or lBournaud etal] (|2004 . 



6. Conclusions 

The present work constitutes the first systematic attempt to try 
to establish TDG candidates among ULIRGs. Extranuclear star- 
forming regions identified as candidates to TDGs have been 
studied in a sample of low-z Ultraluminous Infrared Galaxies 
on the basis of Integral Field Spectroscopy and high angular res- 
olution HST images. 

As previously found in lower luminosity mergers, the present 
work shows that the presence of TDG candidates in ULIRGs is 
common. In particular we have identified twelve condensations 
in 5 ULIRGs (from an initial sample of 9), expanding the lu- 
minosity range of interacting galaxies and mergers showing this 
phenomenon. 

We have characterised the main physical and kinematic prop- 
erties of these candidates. In particular, their dynamical masses 
are in the 2 x 10** - 1 x 10"^ Mq range with typical sizes of 
~ 750 pc. We found that most of these condensations follow 
the relation between effective radius and velocity dispersion of 
lower (globular clusters) and higher (elliptical) mass systems. 

Starbursts have been identified within each TDG candidate. 
Their Hii - like ionisation, relatively high metallicity (12 + 
log(0/H) ~ 8.6), and Ha equivalent widths (34 - 257 A) ai-e 
characteristic of young bursts of star formation (5-8 Myr) with 
stellar masses between 2 x 10'' and 7 x 10*^ Mq. 

Using structural, physical and kinematical information, we 
have discussed the likelihood of survival for these candidates to 
internal as well as to the external forces on the basis of different 
dynamical tracers. Five out of the twelve initial regions present 
a High or High-Medium likelihood of survival as future TDG. 
We identified R in IRAS 15250+3609 as our best candidate. It 
satisfies all our tests, occupies in the r-cr plane a position close 
to dwarf elliptical galaxies, has relatively high relative velocity 
to the parent galaxy and have been found in a system in an ad- 
vanced stage of the merging process. 

Finally, we have outlined the working lines that need to be 
followed in the confirmation of these candidates as well as when 
looking for new ones. That is: IFS observation with high spa- 
tial resolution as those provided by Adaptive Optics assisted IFS 
systems and the systematic search of external regions of star for- 
mation in ULIRGs in an advanced state of the merging process. 
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Fig. 1. Sample of ULIRGs used in this study. First column displays the WFPC2 image in the F814W filter (i.e. HST I-band) for 
reference. Images have been smoothed with a O'.'2-sigma Gaussian and contours have been over-plotted to better see the regions 
under study. Second and third columns present the maps created from the INTEGRAL data for a line-free continuum adjacent to 
Hot and for this line respectively. All images are presented in logarithmic scale. North is up and east is to the left. Regions under 
study have been marked in both the WFPC2 and the INTEGRAL Ha images to better see the matching between images from both 
instruments. 
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Fig. 2. Representative observed spectra of the main extranuclear line-emitting regions. Note that the horizontal axis does not repre- 
sent a continuous range in wavelength, but two subsets of the full INTEGRAL spectrum corresponding to H/6-[Om]AA4959,5007 
and [O i]/l6300 - Hc!'+[N ii] AA6548,6586 - [S ii]/l/l6717,6731 spectral regions, respectively. Flux units are arbitrary and for the shake 
of clarity, offsets have been appUed to the spectra. 
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Fig. 3. Position of the different regions in the diagnostic diagrams proposed bv lVeilleux & Osterbrockl (Il987h . The solid lines show 
the empirical border s found by these aut hors between the different ionisation types while the dashed lines show the theoretical 
borders proposed by iKewley et alJ (l200lh to delimit the area where the line ratios can be explained by star formation. The used 
colour/symbol code is the following: IRAS 08572+3915, cyan/cross; IRAS 12112+0305, green/circle; IRAS 14348-1447, ma- 
genta/square; IRAS 15250+3609, blue/triangle; IRAS 16007+37 43, red/asterisk. Typical errors are ~ 0.2 dex. The measured ratios 
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